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ABSTRACT 

Aims. The recently discovered eclipsing binary system TYC 2675-663-1 is a X-ray source, and shows properties in the optical that are similar 
to the W UMa systems, but are somewhat unusual compared to what is seen in other contact binary systems. The goal of this work is to 
characterize its properties and investigate its nature by means of detailed photometric and spectroscopic observations. 

Methods. We have performed extensive V-band photometric measurements with the INTEGRAL satellite along with ground-based multi-band 
photometric observations, as well as high-resolution spectroscopic monitoring from which we have measured the radial velocities of 
the components. These data have been analysed to determine the stellar properties, including the absolute masses and radii. Additional 
low-resolution spectroscopy was obtained to investigate spectral features. 

Results. From the measured eclipse timings we determine an orbital period for the binary of P = 0.4223576 ± 0.0000009 days. The light-curve 
and spectroscopic analyses reveal the observations to be well represented by a model of an overcontact system composed of main-sequence F5 
and G7 stars (temperature difference of nearly 1000 K), with the possible presence of a third star. Low-resolution optical spectroscopy reveals 
a complex H„ emission, and other features that are not yet understood. The unusually large mass ratio of q = 0.81 + 0.05 places it in the rare 
"H" (high mass ratio) subclass of the W UMa systems, which are presumably on their way to coalescence. 



Key words. Binaries: close - Stars: fundamental parameters (classification, colors, luminosities, masses, radii, temperatures, etc.) - X-rays: 
stars 



1. Introduction 

The star TYC 2675-663-1 (Tycho-2 catalog designation; 
IPfog et all EOOOT) . with coordinates a = 20 h 09 m 11.2 s , 6 = 
+32°33'53" (J2000), was observed during 2002-2005 in the 
Johnson V filter with the Optical Monitoring Camera (O MC) 



on board the INTEGRAL satellite dMas-Hesse et all 120031) . as 



a part of a serendipitous program to monitor optical counter- 
parts of ROSAT sources as potential variable objects. By cross- 
correlating the list of variable sources detected with the OMC 
during the first few months of operation with the ROSAT cat- 
alogs, five were found that were potential optical counterparts 
of X-ray sour ces (i.e., that wer e inside the ROSAT 3cr positio n 
error radius; Caballerol 2004 : Caballero-Garcfa et alJ 2006). 
Three of them showed optical variability with modulations 
characteristic of binary systems (periodic variations, with reg- 
ular periods and amplitudes larger than 0.1 mag). TYC 2675- 
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663-1 showed the most striking variability pattern, with very 
irregular and variable behaviour. 



The 
IOMC 



source, which has received the designation 
267500 0078 in the OMC Input Catalogue 
2003) displays a colour index of B— V ~ 0.7 



Domingo et al., _ _ r 

(H0 getall 120001) . corresponding roughly to spectral type 



G, and shows variability typical of a close eclipsing binary 
system. Little further information is available for this star, 
aside from entries in various astrometric catalogs, so we 
undertook photometric and spectroscopic observations in order 
to ascertain the nature of the object and identify the origin of 
its X-ray emission. 

In this paper we present a detailed study of TYC 2675-663- 
1. The photometric and spectroscopic observations are pre- 
sented in Section [2] with a description of the data reduction 
procedures. Section [3] contains the analysis of these observa- 
tions, including the determination of the ephemeris, the stellar 
parameters and distance to the system, and a discussion of the 
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peculiar H ff emission we detect. Finally, in Section |4] we dis- 
cuss the results and the nature of the system. 

2. Observations and data reduction 

2.1. Photometry 

2.1 .1 . INTEGRAL/OMC photometry 

The Optical Monitoring Camera (OMC) is a 50 mm aper- 
ture refractor telescope, co-aligned with the high-energy instru- 
ments on board the ESA INTEGRAL gamma-ray observatory 
dMas-Hesse et al. , 2003). The OMC provides continuous moni- 
toring of up to 100 sources in the Johnson V band over its 5°x5° 
field of view. Telemetry constraints do not permit downloading 
of the entire OMC image. For this reason, windows were se- 
lected around the proposed X-ray/y-ray targets as well as other 
targets of interest in the same field of view. The position of 
these windows is computed automatical ly, based on the source s 
compiled in the OMC Input Catalogue (Domin go et alj,l2003l) . 
which contains around 500 000 targets selected for being po- 
tentially variable in the optical. Only sub-windows of the 
CCD containing those objects, with a size of 11 x 11 pix- 
els (3.2' x 3.2'), are transmitted to the ground. TYC 2675- 
663-1 is located at only 14" from the ROSAT source 1RXS 
J200912.0+323344 (conta ined in the ROSAT All-Sky Survey 
Bright Source Catalogue, IVoges et al. dl999h . with a catalog 
lcr position error of 8"), and this allowed us to monitor it as a 
potential optical counterpart. 

At each INTEGRAL pointing the OMC collects a set of 
images with several different integration times, typically in the 
range of 10 s to 200 s (currently 10, 50 and 200 s), for the pur- 
pose of optimizing the dynamic range and also to minimize 
noise and cosmic-ray effects. For the analysis of TYC 2675- 
663-1 we have used only the images with exposures of 100 s 
and 200 s, in order to secure an adequate signal-to-noise ratio 
for the observations. The brightness measurements used here 
were obtained with the Off-line Scientific Analysis softwarJ^. 
We used a photometric aperture of 5 x 5 pixels, since contam- 
ination by nearby sources was not significant (see Figure Q}. 
The complete folded light curve is shown in Figure [2] 

2.1 .2. Photometry from the Centro de Astrobiologia 
(CAB) 0.5 m telescope 

Ground-based photometric observations were collected with 
the 0.5 m Giordano Bruno robotic telescope of the Centro de 
Astrobiologia (hereafter CAB), located at Calar Alto, Spain. 
This telescope is equipped with a Finger Lakes Instrumentation 
IMG1024S 1024x1024 back-illuminated CCD, and Johnson 
B and V filters. The camera has an image scale of 0.97 arc- 
sec pixeL 1 , which results in a field of view of 16' x 16'. The ex- 
posure times were 30 s and 20 s for the B and V filters, respec- 
tively. All images were corrected for bias current and flat field. 
We derived differential photometry with respect to three com- 
parison stars in the V filter (USNO Bl 1224-0557970, 1224- 
0558767, and 1223-0560307), and an additional fourth corn- 
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parison star in B (USNO Bl 1224-0558836). This allowed us 
to obtain differential photometry for the source with errors less 
than 0.02 mag for both B and V. 

2.1 .3. Photometry from the Observatorio Astronomico 
Nacional (OAN) 1 .52 m telescope 

In order to obtain better photometric precision we carried 
out additional observations with the 1.52 m telescope at the 
Observatorio Astronomico Nacional, also located at Calar 
Alto (hereafter OAN). This telescope was equipped with a 
Photometries Series 200 back-illuminated CCD, and Johnson 
B, V, and / filters. The camera has a scale of 0.4 arc sec pixel -1 , 
with a field of view of 6.9' x 6.9'. The exposure times were 90, 
60 and 30 s for the B, V, and / filters, respectively. We applied 
standard bias and flat field corrections. Differential photome- 
try was derived with respect to three comparison stars, selected 
from objects in the field with B-V indices similar to the target. 
These stars were USNO-B1 1225-0549546, 1225-0549756, 
and 1225-0549576, with magnitudes of B = 14.1, 15.5, 19.3, 
V = 13.45, 14.65, 18.0, and / = 12.7, 13.7, 16.8, respectively. 
We used a photometric aperture radius of 7" for the target. The 
closest star to the target is located at 11" (USNO Bl 1225- 
0549559) and is fainter, thus its contribution to the flux of the 
target is expected to be negligible. 

2.2. Spectroscopy 

2.2.1. High-resolution spectroscopy 

TYC 2675-663-1 was observed for this project with the CfA 
Digital Speedometer ( LathamL 1992 ) on the 1.5 m Tillinghast 
reflector at the Fred L. Whipple Observatory on Mount 
Hopkins, Arizona (USA). This echelle spectrograph coupled 
with its intensified photon-counting Reticon detector delivers 
a single echelle order 45 A wide centred at a wavelength near 
5187 A, and a resolving power of R « 35,000. The main fea- 
tures in this spectral window are the lines of the Mg I b triplet. 
A total of 18 spectra were gathered between 23 Sep 2004 and 5 
Nov 2006, with signal-to-noise ratios ranging from 14 to 30 per 
resolution element of 8.5 km s . ThAr exposures were taken 
before and after each stellar exposure for wavelength calibra- 
tion. The velocity zero-point was monitored by means of sky 
exposures taken at dusk and dawn. 

2.2.2. Low-resolution spectroscopy 

For a more detailed spectral coverage and to better study the 
spectral properties of the object, a 3-hour sequence of obser- 
vations of TYC 2675-663-1 was gathered also with the 3.5 m 
telescope at Calar Alto using the double-beam spectrograph 
TWIN, on 25 Jul 2007 from 22:02:06 (UTC) to 26 Jul 01 :02:24 
(UTC). The wavelength coverage is -4400 A to 5500 A and 
~6000 A to 7100 A for the blue and red channels, respectively. 
The spectral resolution of these observations is 0.54 A pix~', 
or ~25 and 32 km s _1 in the blue and red channels, respec- 
tively. The exposure time for each spectrum was 600 s (reso- 
lution of 50 phase bins per orbital cycle). The median seeing 
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Fig.l. Field of view of the INTEGRAL/OMC target TYC 2675-663-1 (red point) and the ROSAT source 1RXS 
J200912. 0+323344 (blue circle) at its nominal catalog position. 



Table 1. Log of the observations used in this work. 



Time interval 



Type of observation 



Telescope 



16 Nov. 2002 — 1 May 2005 
21 Jul. 2004 — 1 Sep. 2004 
19 Aug. 2006 — 24 Aug. 2006 
23 Sep. 2004 — 5 Nov. 2006 
25 Jul. 2007 



Photometiy 

Photometry 

Photometry 
High-resolution spectroscopy 
Low-resolution spectroscopy 



INTEGRAL/OMC 
0.5 m CAB 
1.52 m OAN 
1.5 m Whipple Observatory 
3.5 m TWIN 




Phase 



Fig. 2. Complete INTEGRAL/OMC light curve folded with the period of P = 0.4223576 + 0.0000009 d and reference epoch 
Tq = 2 453 080.0249 (time of inferior conjunction of the secondary). 



was about 2". Dome flat fields and bias images were taken at 
the beginning of the night, and comparison HeAr lamp spectra 
were taken regularly for wavelength calibration. 

3. Data analysis 

3.1. Ephemeris 

An initial value of the period of TYC 2675-663-1 was derived 
using a method based on the Phase Dispersion Minimization 
algorithm of Stellingwerf (1978). In order to refine this pe- 



riod, we determined individual times of minimum from all our 
photometry data (INTEGRAL/OMC, CAB 0.5 m, and OAN 
1 .52 m) by fitting each eclipse with a fourth degree polynomial. 
A total of 59 timings were obtained (21 primary eclipses and 
38 secondary eclipses), which are listed in Table|2] These were 
then used to establish the final period and reference epoch by 
solving for a linear ephemeris using standard weighted least- 
squares techniques. Primary and secondary minima were ad- 
justed simultaneously, and the orbit was assumed to be circular. 
Given the asymmetry of the eclipses, realistic uncertainties for 
the individual timings are a bit difficult to determine. Instead, 
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Table 2. Measured times of mid-eclipse for TYC 2675-663-1, 
listed with their errors (cr), cycle number (£), and — C resid- 
uals from the adopted ephemeris. 



Table 2. (Continued.) 



HJD— /4UU UUU 


cr 


O—C 


E 


Year 


lelescope 


52595. 58979 


0.01 10 


i a aaa 1 a 

+0.00910 


-1 147.0 


2002.8764 


OMC 


52595.76130 


0.01 10 


-0.03057 


-1 146.5 


O AA1 OTZTO 

2002.8768 


OMC 


CO C AC CO 

52595.99658 


0.01 10 


A f\f\£ A 1 

-0.00647 


-1 146.0 


2002.8775 


OMC 


52596.20278 


0.01 10 


-0.01 145 


-1145.5 


2002.8780 


OMC 


52596.43136 


0.01 10 


+0.00596 


-1145.0 


2002.8787 


OMC 


CO^ 1 O ylOA^IO 

52612.48048 


0.01 10 


i A f\C\C A A 

+0.00549 


- 1 107.0 




OMC 


52617.96213 


0.01 10 


A A AO C T 

-0.00351 


— 1094.0 


AAAA A*>^7^ 

2002.9376 


OMC 


52619.03278 


0.01 10 


i A A 1 1 O A 

+0.01 124 


-1091.5 


AAA1 A A AC 

2002.9405 


OMC 


52620.29992 


0.01 10 


+0.01 131 


-1088.5 


2002.9440 


OMC 


52620.92656 


0.01 10 


i A A A A A T 

+0.00441 


— 1087.0 


AAAA A A C^7 

2002.9457 


OMC 


52621.76625 


0.01 10 


-0.00061 


-1085.0 


AAA1 A A OA 

2002.9480 


OMC 


co £ o 1 on jirt 

52621.97140 


0.01 10 


A AA£ £ A 

-0.00664 


— 1084.5 


OA AO A/IO^ 

2002.9486 


OMC 


52622.18812 


0.01 10 


A AA 1 1 A 

-0.00110 


-1084.0 


O AAO A A AO 

2002.9492 


OMC 


52622.81 180 


0.01 10 


-0.01095 


-1082.5 


2002.9509 


OMC 


rT/"A') A o Art -i 

52623.03091 


0.01 10 


-0.00302 


— 1082.0 


AAAA AC T C 

2002.9515 


OMC 


52623.46205 


0.01 10 


i A AAC^7^1 

+0.00576 


-1081.0 


2002.9527 


OMC 


CO^O.1 OAAAA 

52624.30990 


0.01 10 


i A A A O O A 

+0.00889 


— 1079.0 


aaaa ArrA 

2002.9550 


OMC 


CO^O A CO~7A^ 

52624.52707 


0.01 10 


i A A 1 A on 

+0.01489 


-1078.5 


2002.9556 


OMC 


52624.71628 


0.01 10 


-0.00708 


-1078.0 


2002.9561 


OMC 


52624.92494 


0.01 10 


A AAA^A 

-0.00960 


— 1077.5 


AAAA 

2002.9567 


OMC 


52737.71 162 


0.01 10 


+0.00760 


-810.5 


O A A O ^l/rr 

2003.2655 


OMC 


52746.76104 


0.01 10 


-0.02367 


-789.0 


AAAA AAAA 

2003.2902 


OMC 


52797.69304 


0.01 10 


i A A 1 A O A 

+0.01424 


-668.5 


AAAO (Am 

2003.4297 


OMC 


52798.32169 


0.01 10 


+0.00935 


-667.0 


2003.4314 


OMC 


CO^AO T1A1 c 

52798.73015 


0.01 10 


A A A A C C 

-0.00455 


-666.0 


AAA1 .lOOC 

2003.4325 


OMC 


CO~7AO AC O /I O 

52798.95242 


0.01 10 


i A f\f\£. C C 

+0.00655 


-665.5 


AAAO l^Al 

2003.4331 


OMC 


CO^AA T^A^O 

52799. 16068 


0.01 10 


i A AAO ^ 1 

+0.00363 


-665.0 


AAAA ^AA'l 

2003.4337 


OMC 


52800.62850 


0.01 10 


-0.0068 1 


-661.5 


2003.4377 


OMC 


COOAA O/^AJ 

52800.86694 


0.01 10 


i A AO A A H 

+0.02046 


-661.0 


AAAO A "tO A 

2003.4384 


OMC 


r/iOAI O £ A A H 

52801.26446 


0.01 10 


A A A A ■> O 

-0.00438 


-660.0 


2003.4395 


OMC 


57806 33100 


0.01 10 


-0 006 1 3 


-648.0 


7003 4533 


OMC 


53208.40495 


0.0124 


-0.01662 


304.0 


2004.5542 


0.5m 


53209.46826 


0.0124 


-0.00920 


306.5 


2004.5571 


0.5m 


53220.45425 


0.0124 


-0.00451 


332.5 


2004.5871 


0.5m 


53224.49223 


0.0124 


+0.02107 


342.0 


2004.5982 


0.5m 


53231.42901 


0.0124 


-0.01105 


358.5 


2004.6172 


0.5m 


53231.65903 


0.0124 


+0.00779 


359.0 


2004.6178 


0.5m 


53233.53572 


0.0124 


-0.01613 


363.5 


2004.6230 


0.5m 


53235.44962 


0.0124 


-0.00284 


368.0 


2004.6282 


0.5m 


53250.43605 


0.0124 


-0.01010 


403.5 


2004.6692 


0.5m 


53324.15075 


0.0110 


+0.00320 


578.0 


2004.8710 


OMC 


53324.59421 


0.0110 


+0.02430 


579.0 


2004.8723 


OMC 


53330.06528 


0.0110 


+0.00472 


592.0 


2004.8872 


OMC 


53330.49275 


0.0110 


+0.00983 


593.0 


2004.8884 


OMC 


53336.39682 


0.0110 


+0.00090 


607.0 


2004.9046 


OMC 


53336.82269 


0.0110 


+0.00441 


608.0 


2004.9057 


OMC 


53342.29115 


0.0110 


-0.01778 


621.0 


2004.9207 


OMC 


53342.74087 


0.0110 


+0.00958 


622.0 


2004.9219 


OMC 


53357.09190 


0.0110 


+0.00046 


656.0 


2004.9612 


OMC 



we assigned reasonable initial uncertainties by telescope, and 
then adjusted them by iterations so as to achieve reduced x 1 
values near unity separately for the minima from each data set. 
In this way we established realistic timing errors of 0.01 1 days, 
0.012 days, and 0.006 days for the INTEGRAL/OMC and the 
0.5 m and 1.52 m telescopes, respectively. The resulting period 
and epoch are given by 



P = 0.4223576 + 0.0000009 days 



T = 2 453 080.0249 + 0.0014 (HJD) 



(1) 



(2) 



where the reference epoch To was chosen to be close to the 
mean value of all the timings. We detect no sign of period 
changes over the 3.7-year interval, and we adopt this ephemeris 
for the remainder of the paper. 



HJD-2 400 000 


(T 


O-C 


E 


Year 


Tel 


53966.55612 


0.0065 


+0.00265 


2099.0 


2006.6299 


1.52m 


53967.39777 


0.0065 


-0.00041 


2101.0 


2006.6322 


1.52m 


53967.61149 


0.0065 


+0.00213 


2101.5 


2006.6327 


1.52m 


53968.44541 


0.0065 


-0.00866 


2103.5 


2006.6350 


1.52m 


53969.51250 


0.0065 


+0.00253 


2106.0 


2006.6380 


1.52m 


53970.36807 


0.0065 


+0.01339 


2108.0 


2006.6403 


1.52m 


53970.56628 


0.0065 


+0.00042 


2108.5 


2006.6408 


1.52m 


53971.40245 


0.0065 


-0.00813 


2110.5 


2006.6431 


1.52m 


53971.61272 


0.0065 


-0.00904 


2111.0 


2006.6437 


1.52m 


53972.46836 


0.0065 


+0.00189 


2113.0 


2006.6460 


1.52m 



3.2. Radial velocities and spectroscopic orbit 

Radial velocities from our high-resolution spectra were de- 
termined using the two-dimensi o nal cr oss-correlation algo- 
rithm TODCOR dZucker & Mazehl 1 1994 . with templates for 
each star selected from a library of synthetic spectra based on 
Kurucz (1992) model atmospheres. Template parameters (ef- 
fective temperature, surface gravity, rotational velocity) were 
chosen to match the properties determined for the components 
in the following sections, and solar metallicity was assumed 
throughout. Although the typical precision for the radial ve- 
locities of single sharp-lined stars with this instrumentation is 
about 0.5 km s _1 , in this case the performance is significantly 
degraded because of the very large broadening of the spectral 
lines (which we attribute to rotation) coupled with the narrow 
wavelength range (single echelle order), and the double-lined 
nature of the object, which combined introduce considerable 
line blending. We estimate individual errors around 25 km s _1 
for our measurements of the velocities of the primary and sec- 
ondary components, which are however much smaller than the 
large variations we detect. These measurements are listed in 
Table [3] A circular Keplerian orbit was fit to these velocities 
holding the period and epoch of primary eclipse fixed at the 
values determined earlier. The fit does not include proximity 
effects, given that the observations were obtained mostly at 
the quadratures where those effects are expected to be small 
(especially when compared to the observational errors). We 
adopt these elements for computing the absolute dimensions 
of the components. The resulting orbital parameters are given 
in Table [4] and the observations along with our best fit are de- 
picted in Figure [3] 



3.3. Light curve fitting 

In order to characterize the system and determine its phys- 
ical properties, we have analysed the light curves taking 
into account the information obtained from the spectro- 
scopic solution just described. Modeling of the photometry 
was performed with the PHOEBE package (IPrsa & Zw itter. 
2005 ) (version 29c) , based on the Wilson-Devinney model 
(Wil son & Devinnevl 1 1 97 lb . The main adjustable parameters 
in this model are typically the inclination of the orbit (J), the 
(pseudo-)potentials (Clu Q2), the luminosity of the primary 
(Li) in each passband, the temperature of the secondary (7^2), 
and a phase shift. The light curves alone do not provide any 
constraint on the absolute temperatures, but only on the tem- 
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MJD 52595-52596 MJD 52612-52620 

11.3 i | i | I | I | I | I I 11.3 I i | i | i | i | r | i 




Phase Phase 

Fig. 4. INTEGRAL/OMC light curve of TYC 2675-663-1 in the Johnson V band, folded with the adopted ephemeris. Separate 
time intervals are shown to illustrate the changes in shape due to activity (from upper-left to lower-right). The superimposed 
curves (red) represent a fit without accounting for spots (see text), shown here only for reference. 



perature ratio through the depth of the eclipses. Therefore, the 
primary value T^i was held fixed. We estimated the tempera- 
ture of the primary star from the combined (V - K) colour and 
estimates of the magnitude difference between the components 

2 We adopt V = 11.40 ± 0.05 and K = 10.116 ± 0.023, where V 
comes from the mean INTEGRAL/OMC standard magnitude out of 
eclipse (similar to the Johnson V magnitude). The corresponding un- 
certainty is a conservative estimate based on the short period variations 



derived from the light curves, in an iterative way. The result is 



of the light curve over timescales of a few weeks (less affected by stel- 
lar activity, although those effects are still present, as seen in Figure|4j- 
The K magnitude is taken from the 2MASS catalog i Cutri et all 
2003). We adopt also reddening correctio ns of E(V- K) = 0.16 and 
E(B-V) = 0.065, as well as A(V) = 0.2 dFitzpatrickLll999[V The tem- 
peratu re is based on the (V - K)-T c g relationship from lMasana et al] 
(2006). 



6 



M. D. Caballero-Garcia et al.: The new binarv system TYC 2675-663-1 

_, . . . , . , n V band 



200 



E100 



£100 



-200 
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Fig. 3. CfA radial velocity measurements of TYC 2675-663- 
1 along with our best-fit model, folded with the period P = 
0.4223576 + 0.0000009 d and T Q (HJD) = 2453 080.0249 + 
0.0014 (corresponding to phase 0.0). 



Table 3. Radial velocities obtained for the primary (1) and the 
secondary components (2). 



HJD-2400000 


RVi (km r J ) 


«V 2 (km s _1 ) 


53271.6872 


+ 180.87 


-125.13 


53272.7286 


-103.84 


+ 177.81 


53273.7860 


+ 188.05 


-158.08 


53275.6659 


-141.83 


+ 190.85 


53278.6415 


-137.75 


+208.62 


53280.7467 


-162.68 


+203.43 


53281.6124 


-121.20 


+235.58 


53281.7662 


+ 156.63 


-119.48 


53282.6452 


+ 182.05 


-159.53 


53301.6292 


+ 159.51 


-146.73 


53308.6291 


-71.43 


+ 160.49 


53684.7073 


+201.79 


-176.26 


53686.6018 


-107.92 


+226.31 


53690.6289 


+ 188.68 


-152.36 


53691.6825 


-75.28 


+212.21 


53693.5844 


+ 161.95 


-208.17 


53873.9388 


+ 157.45 


-197.33 


54044.5790 


+ 191.78 


-159.58 



Table 4. Spectroscopic orbital parameters of TYC 2675-663-1. 



Parameter 


Value 


q = M 2 /Mi 


0.81 ±0.05 


a sin i (R e ) 


2.85 ± 0.08 


v y (km S" 1 ) 


27 + 4 


e 


°1 


K x (kms- 1 ) 


153 + 6 


K 2 (km s -1 ) 


189 ±7 


Mi sin 3 i (M ) 


0.97 ± 0.09 


M 2 sin 3 i (M ) 


0.78 ± 0.07 



" Circular orbit assumed in the model. 




B band 




I band 




Fig. 5. Light curves of TYC 2675-663-1 obtained with the 
OAN telescope in the V, B, and / bands (top to bottom), folded 
according to the adopted ephemeris. This photometry is strictly 
differential (magnitudes shown have an arbitrary zero point). 
The curves represent our best-fit solution for an overcontact 
configuration, including the effect of spots. 



T e ff,\ = 6480 + 180 K, where the error is a conservative esti- 
mate based on the scatter of values found from the colour in the 
iterative process of fitting the OAN light curves. 

Limb -darkening coeffici ents were interpolated from the ta- 
bles by Ivan Hammd dl993l) . using the square-root law. The 
gravity brightening exponents (0.3 for both stars) and bolomet- 



ric albedos (0.5) were fixe d at the value s appropriate for stars 
with convective envelopes (ILucvll 19681) . The mass ratio q and 
projected semimajor axis a sin i were adopted from the spec- 
troscopy. Because of the higher quality of the BVI photometry 
from the 1.52 m OAN compared to that from the 0.5 m CAB, 
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Table 5. Light curve solutions based on the OAN photometry, 
for an overcontact and semi-detached configuration. 



Parameter 


Overcontact (adopted) 


Semi-detached 


in 


65.2 + f 


62 90 +2 211 


T-eff.2 (K) 


5543+™ 


5724+jf 


loggi 




4 49^+0.015 


l°gg2 


4 36 +0 06 


4 40 +0 04 


ill 


Q Q7+0.06 

-"•-"-0.05 


3.43 ± 0.08 






T QQ7+0.018 

J - Jyl -0.030 


Li/(Li +L 2 + L 3 ) B 


0.71 ±0.09 


0.685 ± 0.006 


L,/(Li +L 2 + L 3 ) V 


0.680 ± 0.029 


0.66 ± 0.05 


UKU +L 2 + L 3 ), 


0.6272 ± 0.0016 


0.613 ±0.016 




0.035 ± 0.018 


0.0263 ± 0.0027 




0.0399 ± 0.0006 


0.031 ±0.009 


hi 


0.0489 ±0.0011 


0.037 ± 0.003 


>> 

Ppolc 1 _ 


0.38 ± +0.05 


0.37 ± 0.05 


Ppoint 1 


- 


0.52 ± 0.05 


Pside 1 


0.40 ± 0.05 


0.39 ± 0.05 


Pback 1 


0.44 ± 0.05 


0.42 ± 0.05 


Ri (Rot 


1 99+0.05 

1,z - -0.03 


1.190 ±0.005 


Ppolc 2 


0.35 ± 0.05 


0.34 ± 0.05 


Ppoint 2 






Pside 2 


0.36 ± 0.05 


0.36 ± 0.05 


Pback 2 


40 + n 

U.H-U It U.UJ 


U.J7 It U.UJ 


Rl (Re) 


1 11+0.05 
1 - 11 -0.03 


1 10 +002 
i - iu -o.oi 


^boii (mag) 


o 01 +0.18 
J - S1 -021 


3.86 ±0.12 


M bo i2 (mag) 


4 69+° 08 
Spot parameters 


4 57+0.04 
^ ,J '-0.16 


Colatitude 2 (deg) 


20 


20 


Longitude 2 (deg) 


120.0!* 5 


135.2 + ' 9 2 


Radius2 (deg) 


57 +47 ' 


74+12 

' -23 


T'spot 2/^surface 2 


0.80 


0.80 


^Acr 2 (V)/N 


0.001136 


0.001240 


VsAo- 2 (e)/yv 


0.000910 


0.001006 


VSAo- 2 (7)/JV 


0.002939 


0.003108 



" Third-light values at phase 0.25. 

b Relative radii in different directions, in units of the semi-major 
axis. 

c Volume radius, defined as the radius of a sphere with the same 
volume as the limiting Roche lobe. 



Fig. 6. Configuration of the system at different phases, showing 
the location of the spot on the secondary. From top to bottom, 
the phases are approximately 0.00, 0.25, and 0.75. The size and 
separation of the stars are rendered to scale. 



and the relatively short time span of the OAN observations that 
makes them less vulnerable than the INTEGRAL/OMC pho- 
tometry to variability in the light curve, as described below, 
we use only the OAN data in the following to determine the 
light curve parameters. The three passbands were fitted simul- 
taneously. The INTEGRAL/OMC measurements were used as 
a consistency check. 

It is obvious from the light curves (Figure [4]i that the shape 
of the modulation changes continuously with phase, as in the 
classical W UMa systems, with no clear beginning or end- 
ing of the eclipses. This strongly suggests significant defor- 



mation of the components, and perhaps some degree of con- 
tact. Consequently, we performed fits in both the semidetached 
mode (mode 4 in the Wilson-Devinney nomenclature, with the 
primary filling its Roche lobe), and also in the overcontact 
mode (mode 3), which is appropriate for W UMa-type sys- 
tems and overcontact systems that are not in thermal contact. 
In the former case the potential of the primary is set to the 
value of its Roche lobe potential for the adopted m ass ratio 
(q = 0.81 + 0.05): Oi = 3.43 + 0.08 dWilsonlll99 A In the 
overcontact mode the potentials are constrained to be the same, 
but we do not force there to be in thermal contact, in other 
words, the surface brightness of the two stars can be different 
even though they might be in geometrical contact. 



3 The corresponding secondary Roche lobe potential is Q.2 
0.14, assuming synchronous rotation. 
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Table 6. Light curve solutions without accounting for spot ef- 
fects, based on the INTEGRAL/OMC and OAN photometry. 
These results are for comparison purposes only, and are not 
our final adopted values. 



Parameter 


Fits to INTEGRAL/OMC 


Fits to OAN 




data 


data 




68 -" 


64 8 +L , 

" -1.4 


* en, i V 1 -*-/ 


5400 + ?^ 


5350+,,?, 

JJJ "-140 


log j?i 

fa o 1 


4.32+0-i^ 

-0.16 


4.36 ± 0.09 


log go 


4.30+ ,!' 

-0.15 


4.34 ± 0.09 




o or +0.07 
-0.09 


t ^9+0.03 
-0.02 






" I 


1 ,1 1(1 A + /,i 4- Li\d 




0.737 ± 0.024 


Li /(Lt + L-> 4- Li ^ i/ 


0.70 ± 0.09 


0.7172 ± 0.0020 


Li /(L\ + Li + Li)t 




0.664 ± 0.012 






0.031 ±0.006 




061 +008 

u - uu -0.015 


0.0433 ± 0.0010 






0.0483 ±0.010 


b 

Ppolc 1 _ 


0.39 ± 0.05 


0.39 ± 0.05 


Ppoint 1 






Pside 1 


0.42 ± 0.05 


0.41 ± 0.05 


Pback 1 


0.46 ± 0.05 


0.45 ± 0.05 


Rl (Ret 


1 26 +0 - 06 


1.25 ± 0.05 


Ppolc 2 


0.36 ± 0.05 


0.35 ± 0.05 


Ppoint 2 






Pside 2 


0.38 ± 0.05 


0.37 ± 0.05 


Pback 2 


0.42 ± 0.05 


0.42 ± 0.05 


Rl (Re) 


1 16 +005 

1-1 -0.08 


1.15 ±0.05 


A/boii (mag) 


o 74+0.25 
J - ' -0.24 


3.76 ±0.21 


M bo[2 (mag) 


471 +0.7 
H - /J --0.25 


4 71+0.27 
^- ' 1 -0.21 



" Third-light values at phase 0.25. 

* Relative radii in different directions, in units of the semi-major 
axis. 

Volume radius, defined as the radius of a sphere with the same 
volume as the limiting Roche lobe. 

Table 7. Absolute dimensions of TYC 2675-663-1. 



Parameter 



M(M ) 
R(Ro) 
T eB (K) 

logs 
L(L ) 
M bol (mag) 
BC V (mag) 
My (mag) 



Primary 

1 9Q+0.11 

i 99+0.05 
^ -0.03 

6480 ± 180 

4 38 +0 05 

-0.032!°;™ 
3-84!gj° 



Distance (pc) 



360 



Secondary 

1 - <JH -0.04 

i ir 005 

1 ' 11 -0.03 

5543+5°° 
4 36 +0 ° 6 
1 04 +06 

i ' UH -0.07 

4 59+0.08 

^' u7 -0.50 

_n 14Q+0.080 

u - 1M ^-0.005 

4 S4+0.08 

"T313 

-50 



Initial solutions indicated a secondary temperature T e g t 2 in 
the range of 5100-5700 K, but also gave rather low values of 
the inclination angle 57 deg) as a result of the relatively 
shallow eclipses, leading to absolute masses of Mi « 1 .62 M Q 
and Mi ~ 1.31 M Q . Such large masses, which are typical of 
late A- and mid F-type stars, respectively, would imply tem- 
peratures that are considerably higher than we estimate based 
on the colours (see above). Experiments in which we included 
third light (£3) in our solutions revealed that it is statistically 



different from zero (> 3<x), providing a plausible explanation. 
We obtained (3 values of 2-3% in B, 3-4% in V, and 4-5% 
in /, which suggests a red object contaminating the photome- 
try. These solutions yielded considerably larger inclinations, as 
expected, by about 5-9°. We do not see any clear evidence of 
a third star in our CfA or TWIN spectra, although the spectro- 
scopic material is inadequate for detecting such faint signatures 
so it does not rule them out. In the following we have chosen to 
include third light in our light curve modelling, on the basis that 
it is statistically significant and provides for a more consistent 
overall solution. 

There are abundant signs of chromospheric activity in 
TYC 2675-663-1, manifested in the form of an occasional 
O'Connell effect, irregular shapes of the minima, and occa- 
sional flares (see Figure |4|. There is also significant variabil- 
ity in the shape of the light curves, on typical timescales of a 
few weeks. The O'Connell effect is apparent in the different 
light levels at the quadratures (AV, AB * 0.1 mag), and its 
reality is demonstrated by the fact that it is seen in the light 
curves from the 0.5 m CAB telescope (not shown) during one 
of the same time intervals covered by the INTEGRAL/OMC 
photometry (Figure|4|, thus ruling out instrumental errors. This 
strongly suggests spottedness in one or both stars, which must 
be considered in order to avoid biases in the geometric ele- 
ments. PHOEBE allows the spot effects to be modelled, by 
assuming one or more uniform-temperature circular features 
parametrized by four additional adjustable variables: the longi- 
tude, colatitude, angular size, and temperature contrast factor 
(T spot /T 'star)- However, because of well-known degenerac ies in 



Eker 



fitting for s pot parameters from light curves (see, e.g. 
Il996lll999h . and the limited quality of our observations, it is 
not possible to discern precisely and unambiguously where the 
spots are in this system. For this work we have arbitrarily cho- 
sen to place a single spot on the secondary star, and we have 
chosen its location to be near the pole (colatitude b = 20°), 
on the grounds that similar locations are often seen in other 
active and rapidly-rotating stars studied, e.g., with Doppler 
imaging techniques (see, e.g. Strassmeier 2009 and references 
therein). We have also fixed the temperature contrast factor to 
7^pot/7"stai- = 0.80, a value similar to that used in other studies. 
Experiments with a spot location near the equator produced fits 
of about the same quality, and also gave similar values for all 
geometric and radiative quantities, within the errors. The fits 
converged to a large spot covering ~20% of the surface of the 
secondary star. 

Solutions with third light and spots using the OAN pho- 
tometry were carried out both in the semi-detached mode and 
in the overcontact mode. The fits consistently indicated con- 
tact between the components and a slight preference for the 
overcontact configuration, but gave otherwise very similar re- 
sults for all parameters, the differences being well within the 
errors. As seen in Table|5]the potentials of the two components 
in the semi-detached solution appear nearly identical (within 
the errors), and equal to those from the overcontact solution. 
This fact agrees with the definition of the overcontact mode, in 
which the components have the same potentials. Therefore, we 
adopt in the following the overcontact results. These solutions 
are presented in Table |5j and are shown in Figure [5] together 
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with the OAN photometry. The values of the parameters corre- 
spond to the results from fitting simultaneously the B, V, and / 
light curves, and the errors represent the dispersion of the so- 
lution with respect those from fits of the monochromatic light 
curves. A graphical representation of the configuration of the 
system and the location of the spot is shown in Figure [6] 

As a consistency check, we produced an additional fit using 
the INTEGRAL/OMC data in the adopted overcontact mode. 
However, because of the longer time coverage of these data 
compared to the OAN light curves, and the significant spot 
variability on relatively short timescales, we did not consider 
spots in this solution (i.e., this is not an optimized fit, and not 
our adopted fit). For the purpose of comparison, we fitted the 
same spotless model to the OAN light curves, and the results 
are listed in Table[6] There is very good agreement between the 
two data sets, despite the simplified modelling. This solution 
is over-plotted with the INTEGRAL/OMC data in Figure H A 
comparison for the OAN data between the spotted model in 
Table [5] and the unspotted model in Table [6] shows that the dif- 
ferences are rather small, and thus that the absolute dimensions 
of TYC 2675-663-1 are unlikely to be much affected by the 
uncertainties in the spot modelling. 

3.4. Absolute dimensions 

The combination of the light curve and spectroscopic param- 
eters leads to the physical properties of the system presented 
in Table [7] The temperature difference between the stars is 
~940 K, and the individual values correspond to spectral types 
of approximately F5 and G7 for the primary and secondary, re- 
spectively. Using these individual temperatures, we infer bolo- 
metric corrections BC V of -0.032+° °™ and -0.149+°,-°$? for the 
primary and secondary stars, respectively, based on the tabula- 
tions bv iMasana et al.l ([2006). The errors are propagated from 
the temperature uncertainties. With the bolometric magnitudes 
that follow from the temperatures and radii, we derive absolute 



visual magnitudes My of 3.84+Q2° and 4.84+°°y. The magni- 
tude difference is in good agreement with the value inferred 
from the flux ratio between the components from the fits (i.e., 
My 2 -M V i = 2.51og(Z4/L 2 ) = 0.96!°;^; see Table E). 

The apparent visual magnitudes of the primary and sec- 
ondary follow from the out-of-eclipse brightness and the light 
ratio, and are 1 1.80 + 0.08 and 12.80 ± 0.08, res pectively. With 
a visu al extinction estimate of A(V) = 0.2 from lDrimmel et al. 
d2003h . and the absolute magnitudes given above, we infer a 
distance to the binary system of D 



360!™ pc. 



3.5. X-ray emission 

The sky position of the optical source TYC 2675-663-1 is 14" 
far from the ROSAT source 1RXS J2 00912.0+323344 which 
has a lcr positional uncertainty of 8" dVoges et al. , ll999l) . This 
implies a formal coincidence likelihood between the optical 
and the X-ray locations of < 13%. The reality of this X- 
ray source is confirmed by its detection by the XMM-Newton 
satellite (under the designation XMMSL1 J200 910.0+323358) 
at similar coordinates as the ROSAT source (IFrevberg et al. 



2006) and with a similar positional uncertainty^. The count 
rate and hardness ratio of the source measured by ROSAT are 
0.123 + 0.016 counts s" 1 and 0.82 + 0.08, respectively. This 
translates to an integrated fl ux of 1.56xl0~ 12 erg s~ 1 cm~ 2 in 



the 0.1-2.4keV energy band dSchmitt et al. , 1995 ). Assuming 



a distance to the source of D — 360!™ pc then gives an intrin- 
sic luminosity of Lx = (2.4+_qJ)x10 ergs -1 . Considering the 
bolometric luminosities of the components of Lboi = 2.4!^ L 
and 1 .04!°q 7 L Q for the primary and the secondary, we derive 
an X-ray to total luminosity ratio of log Lx/L\, \ = -2.74!^^, 
which is typical of binary systems w ith rotationally induced 
activity (see, e.g., Messi na et al ] l2003h . 



3.6. Analysis of spectroscopic features 

Figures [7] and [8] display the sequence of 14 normalized low- 
resolution spectra from the TWIN spectrograph in the range 
from 4800 to 4950 A and from 6400 to 6700 A, showing the 
most interesting features. Among the most prominent are the 
strong Balmer lines in emission (H a and H^ at A = 6563 A and 
4861 A, respectively) and the He i line centred at 6680 A. The 
H tl , emission is complex and is discussed below. 

In the red spectra (Figure [8) there is a clear absorption 
feature at around 6490 A. If this has a photospheric origin, it 
is likely to be a blend composed mainly of Ti n, Ca i, Fe i, 
and Ba n lines, which are the ones with the largest oscillator 
strengths in this spectral range. To illustrate this, in Figure [9] 
we show the synthetic spectrum of a star similar to the primary 
component of TYC 2675-663-1 (r eff ~ 6000 K, logg = 4.0, 
v sin i = km s _1 , and an assumed [M/H] = 0.0), together with 
the result of applying a rotational broadening of 144 km s _1 , 
which is our estimate for the primary star assuming spin-orbit 
synchronization. This washes out the narrow lines, producing a 
single broad feature, as observed. The steady character of this 
broad absorption suggests its origin in a third star, perhaps a 
distant companion to the binary system. This may perhaps ex- 
plain the third light contribution inferred from the fits of the 
light curves in Section [331 Nevertheless, more complete spec- 
troscopic observations are needed in order to confirm the likely 
stellar origin of this absorption complex. 

3.6.1. The Hq, emission line 

Very broad Hq, emission (Av w 1200km s ) is detected in 
the spectra of this binary system. Phenomenologically, the Hq, 
complex can be described by the presence of what appear to be 
two broad and asymmetric emission lines centred at 6553 A and 
6568 A, plus a steady absorption dip between them (P-Cygni 
profile). This emission has a complicated structure, with sev- 
eral components that seem to be moving with orbital phase. In 
order to better understand this complex evolution with phase 
we have fitted the profile with additional components. After 



4 The quoted positions in the slew catalog have a 1 sigma error of 
8", which for point sources is dominated by the accuracy of the at- 
titude reconstruction during the slews. The statistical position error 
quoted in the catalog has a mean of ~4" (1 image pixel) for non- 
extended sources. 
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Fig. 7. Blue spectra of TYC 2675-663-1 obtained with the 
TWIN spectrograph (4800-4950 A) showing the H^ emission 
complex centred at Aq = 4861 A. Orbital phases for each spec- 
trum are as indicated in Figure [8] 

trying a number of different combinations, we find that the be- 
haviour of the H„ region of the spectrum can be adequately 
described with five components (see Figures [TOlandfTTIi: one in 
absorption and four in emission. The absorption line (line 5) is 
steady, with a radial velocity of RV ~ -200 km s _1 . 

Similar outflows have been associ ated with stella r winds 
in high mass or pos t-AGB stars (IVarricatt et all |2004; 



Chatv & Filliatre, 2005), with velocities similar to those we 



Smolinski et al .119931) . However, in the 



find (e.g., -105 km s~ 
case of the present binary system TYC 2675-663-1, which is 
composed of two main sequence stars of late spectral type (and 
likely a third star of similar or later spectral type), the absorp- 
tion is unlikely to be due to stellar winds. Its nature is presently 
unclear. 

In Figure QT| two of the Gaussian components (lines 3 and 
4) appear to follow the orbital motion of the primary star, with 
a maximum velocity excursion at quadrature (0 = 0.25). The 
velocities at quadrature are RV « 150 km s and RV ~ 575 
km s _1 , respectively. Given the previously measured orbital 
semiamplitudes of 153+6kms~ 1 and 189+7kms~' for the pri- 
mary and secondary, respectively, we tentatively identify line 3 
with the surface of the primary star. The emission represented 
by line 4 comes from a region that is moving 3—4 times faster. 
Elucidating the origin of these wings and the absorption fea- 
ture of the Hq, complex will most likely require further spec- 
troscopic observations that are beyond the scope of the present 
work, in order to obtain complete phase coverage over a full 
orbital cycle. 

4. Discussion 

Our analysis of the close eclipsing binary TYC 2675-663-1 
has revealed properties similar in many respects to those of 
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Fig. 8. Red spectra of TYC 2675-663-1 obtained with the 
TWIN spectrograph (6400-6700 A) showing the prominent 
Hq, and He i emission features, centred at Aq = 6563 A and 
Ao = 6680 A, respectively. In this figure the wavelength shift 
due to the orbital motion of the components is approximately 
4 A. Orbital phases are labelled. 




X (A) 



Fig. 9. Synthetic spectrum corresponding to a star similar to 
the primary of TYC 2675-663-1, with T eff * 6000 K, logg = 
4.0, [M/H] = 0.0, and no rotational broadening (red line). 
Superimposed is the same spectrum broadened with a rotation 
profile of 144 km s _1 (black line), estimated for the primary in 
TYC 2675-663-1. 



the W UMa systems, which are characterized by having short 
orbital periods (0.2-0.8 d) and an overcontact configuration, 
and are composed of F-K stars sharing a common envelope 
that thermalizes the stars. This leads to near equal depths for 
the eclipses. TYC 2675-663-1 displays eclipses with clearly 
different depths, which would imply non-thermal equilibrium. 
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Fig. 10. Enlarged view of the H a complex from Figure [8] with 
our 5-component fits (one absorption component, and 4 emis- 
sion components; see text). Orbital phases are labelled. 
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Fig. 11. Radial velocity variation of each line component fitted 
to the H„ complex in Figure |T0j as a function of orbital phase. 
The red line represents the absorption component, and the other 
four lines are in emission (see text). 



Nevertheless, the stars have spectral types of F5 and G7 in 
the typical range for the W UMa class, and appear to be 
in geometrical contact despite the temperatures being differ- 
ent by approximately 940 K. The large mass ratio we derive, 
q = M 2 /Mi = 0.81 ± 0.05, identifies the object as an H-type 
W UMa vari able ("H" for high m ass ratio), a subgroup first 
proposed by ICsizmadia et al.l (120041) . In these objects the en- 
ergy transfer rate appears to be less efficient than in other types 
of contact binaries. They also show excess angular momentum, 
which can be understood as resulting from t he first s t ages o f the 
interaction between the components (see iLi et all |2004, and 



references therein). The overall properties of TYC 2675-663- 
1 suggest a system at or near contact and perhaps in the early 
stages of evolution toward a state of full geometrical and ther- 
mal contact (coalescence). 

W UMa systems such as the binary studied here typi- 
cally show increased chromo spheric activity du e to the rapid 
rotation of the components (I Apple gate L Il992[). and this i s 
usual ly accompanied by soft X-ray emission dMessina et al 



2003). This chromospheric activity is manifested by the un- 



equal brightness at quadratures, asymmetrical minima, and er- 
ratic flares shown in the light curves, which are collectively 
referred to as the O'Connell effect. TYC 2675-663-1 shows 
clear evidence of each of these, including the X-ray emis- 
sion. The presence of the He i A 6680 A emission line, as we 
see in the spectra of TYC 2675-663-1, is often interpreted 
as evidence of the interaction between the components of 
the binary sy stem, in the form of winds or streams of mat 



ter (iGrahaml 1 19921: iGreelev et al. 

IScholz & Javawardhanall2006h . The H a region of the spectrum 
is complex, with emission and absorption components that are 
not yet fully understood. 

Our spectroscopic and photometric observations of 
TYC 2675-663-1 have provided a first picture of this close bi- 
nary system in which the two stars differ greatly in tempera- 
ture, but are in a near-contact configuration perhaps leading to 
coalescence. It appears to be a rare example of this class of 
W UMa objects. There is circumstantial evidence for a third 
star in the system. It is hoped that this study will be helpful in 
the development of theories to understand the early behaviour 
and evolution of W UMa systems. 
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